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ABSTRACT 

Constraining the astrophysical nature of Ultra-Luminous X-ray (ULX) sources, which have 
X-ray luminosities exceeding 10 39 erg s -1 , has been elusive due to the optical faintness of any 
counterparts. With high spectral resolution observations in the ~ 10-30/zm wavelength range 
we have conducted an experiment to study six ULX sources in the NGC 4485/90 galaxy pair. 
We have found that five of the six ULXs, based on mid-infrared spectral diagnostics, show the 
characteristic higher ionization features that are found in AGN. The sixth source, ULX-1, is 
consistent with being a supernova remnant. The chief infrared spectral diagnostics used are the 
ratios of [S III] /[Si II] vs [Ne III]/[Ne II]. In two instances fits to the continuum and poly aromatic 
hydrocarbons (PAH) features also indicate higher dust temperatures that are characteristic of 
accreting sources. Overall, however, we find the continuum is dominated by stellar processes, 
and the best diagnostic features are the emission lines. High spectral resolution studies in the 
mid-infrared thus appear to show great promise for determining the astrophysical nature of ULXs. 

Subject headings: galaxies: general — ULRGS: individual(NGC 4485, NGC 4490) 
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1. Introduction 

Ultra-Luminous X-ray sources (ULXs) are gen- 
erally denned as extra-nuclear galactic X-ray 
sources with X-ray luminosities Lx > 10 39 erg s -1 , 
a luminosity that is seldom reached by any Galac- 
tic (Milky Way) X-ray binary (XRB). If ULXs 
emit their X-rays isotropically, by Eddington ar- 
guments, the central object could be a black hole 
(BH) with mass M > 5 M Q , and for Edding- 
ton ratios of ^10~ 2 — 1CU 1 , the implied mass of 
the central object is ~50— 500 M©. Therefore, 
some ULXs could be accreting intermediate-mass 
black holes (IMBHs. IColbert fc Mushotzkvlll999t 
Ptak fc Griffiths! 1 19991 iMiller fc Colbertl l2004 . 
The majority of ULXs, however, are likely ac- 
creting stellar mass black holes (e.g., ~ 1-10 M ) 
as there appears to be a strong association be- 
tween ULX sources and s tar formation in galaxies 
(|Zezas k Fabbianol l2002f) . One popular model 
claims that ULXs are merely stellar-mass BHs 
in high-mass XRBs, and mild beaming of their 
X-ray luminosities produces the observational il- 
lusion of extremely high X-r ay luminosities (e.g., 
King et al.|[200ltlKing ll2006h . Hundreds of ULXs 
are known at the present time (e.g., Roberts k 
Warwick 2000, Colbert k Ptak 2002, and Chan- 
dra ULX catalogs by Ptak et alQ, Swartz et al. 
2004, and Liu & Bregman 2005). 

Accreting binaries such as these could be pre- 
cursors to eventual gravitational radiation events. 
Further, ULXs are important as they often dom- 
inate the to tal observed X-ray flux in starburst 
galaxies (see Colbert et al. 2004 ). Their relation- 
ship with star formation also appears to extend 
to cos mologically interesting (z > 0.1) look-back 



times ( Hornschemeier et al.l l2004t iLehmer et al 



20061 ) and may provide a cosmological tool. 

Unfortunately, all but the very brightest ULX 
sources are too faint for detailed X-ray spec- 
troscopy studies (the best cases must rely upon 
resolving power o f R= 60 at 6.4 k ev, available 
with X-ray CCDs: IMiller et al.ll2004l ). 

There is perhaps hope for ascertaining the na- 
ture of ULX sources by studying their interactions 
with their surroundings, such as gas, dust and 
stars; such information may be captured u si ng in- 
frared spectral diagnostics. iGenzel et al. (1998) 



were the first to show that ionization-sensitive in- 
dices based on mid-infrared ratios were helpful for 
identifyi ng the nature of heavily obscured nuclear 
sources (IGenzel et al. 1998: lLaurent et al.l 2000; 



Sturm et alj2002l|Peeters, Spoon, fc Tielensfl2004bf ) 
We have thus taken advantage of the resolving 
power (R ~ 600) of the Infrared Spectrograph 
(IRS) on board the Spitzer Space Telescope (SST) 
to calculate highly accurate ionization diagnostic 
ratios. 

The nearby galaxy pair NGC 4485 /90 provides 
an e xcellent opportunity for studying ULX sources 
(see iRoberts et ail l2002h . There are six ULXs in 
NGC 4485/90 (D = 7.8 Mpc. lTullvlfl998h . far more 
tha n is typical (< 2 in most cases, and 0.2, on aver - 
age; IColbert fc P"tai3l2002t iPtak fc Colbert |2004 . 
The unprecedented sensitivity and angular resolu- 
tion afforded by the Spitzer Space Telescope's In- 
frared Spectrograph (IRS) coupled with the prox- 
imity of NGC 4485/90 allow us to carry out this 
important experiment. 

2. IRS spectroscopy: Observing strategy 
and Reduction 

We have performed spectral mapping of the sys- 
tem NGC 4485/90 by using the mapping mode of 



the IRS (Ho uck et al J 12004 ) on board the SST 



1 URL: http://www.xassist.org 



Rectangular regions of ll'.'l x22'/3 and 4'/7xll'/3 
(i.e. zones of — 0.42 x 0.84 kpc and ~ 0.18 x 
0.42 kpc at the system) were mapped using both 
the short-high (SH, ~ 9.9 - 19.7 fim) and long- 
high (LH, ~ 19.5 — 38 /mi) slits, respectively Six 
rectangular regions centered at the positions of the 
ULXs (Table [T]) were covered using a 1 x4 and 1x7 
grid array for the long and and short slits of the 
IRS. In addition two positions in NGC 4490 that 
are not near ULXs were observed with the same 
configuration as the one used for the ULXs. The 
integration time was 240 s at each position, and 
the step size is half the slit width perpendicular to 
the slit. This configuration allows net exposures 
of > 480 s for rectangular sizes 16.7x22.3" and 
14.1x11.3" for long and short high slits, respec- 
tively. 

Sets of 8 exposures were taken for each po- 
sition with the IRS. The pre-reduced basic cal- 
ibrated data (bed) files were obtained from the 
Spitzer Science Center (SSC) data base. In order 
to extract the spectra for the SH and LH modes, 
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we used SMART ( Spectroscopy Model ing Anal- 
ysis and Reduction. iHigdon et al. 2004 ) which is 
publicly available software provided by the SSCJE 

3. Infrared Spectral Diagnostics 



We have used line i denti fications from Smith et al. 
( 2004l ). lArmus et al.l(|2006l ) and NISlFI (|Martin et al 
19951 ). We have combined both modes by re- 
scaling the LH mode spectra to the SH mode 
spectra by adding or subtracting a constant and 
then re-binning to a resolution of 0.0287 /tin/pixel. 
In Fig. [T] we show the combined spectra from all 
the positions. 

As we can see from Fig. [IJ some of the line fea- 
tures are present or stronger in the ULX region 
than in the comparison region spectra. The most 
striking examples are features such as [S IV] 10.51, 
[Ne II] 12.81, [Na I]14.6, [Ne III] 15.55, [S III] 18.71, 
and [S III]33.48 fim that seem to be significantly 
stronger in the comparison regions and ULX- 
1 (which is thought to be a SNR rather than an 
accreting binary), than in the rest of the sources. 
In the opposite way, the [Si II] 34.82 /im line is 
stronger in ULX-2 to ULX-6. 

The purpose of this work is to generate some 
emission line diagnostic diagrams similar to 
commonly used optical emission line diagrams, 
which have proven useful in separating accretion- 
dominated sources from those dominated by stel- 
lar sources (e. g.. iBaldwin. Phillips fc Terlevich 



1981c iKauffmann et al.ll2006l . and references therein) 
A popular mid-infrared diagnostic first put forth 
in an ISO spectroscopic surv ey of ultra-luminous 
infrared galaxies (ULIR Gs) by Genzel et al.l dl998l) . 
and lat er explored by IPeeters. Spoon, fc Tielensl 
(|2004bl ). plots the emission line ratio [O IV] 25.91 
/im / [Ne II] 12.81 /im versus the strength of a 
mid-infrared poly aromatic hydrocarbons (PAH) 
features. These two lines however, are rela- 
tively weak and difficult to detect in lower lu- 
minosity systems. Taking advantage of the 
stro ng low ioniz a tion l ine emission of [Si II]34.82 
/xm, Dale et al. (2006) suggested alternative di- 
agrams to distinguish AGNs from SF regions 
such as [Ne III]15.55/im/[Ne II]12.81/im and 
[S III]33.48/mi/[Si II]34.82/im. These ratios 



2 http: / /ssc. spitzor.caltech.edu/archanaly/contributcd/smart/ 
3 Available at http:/ /physics. nist.gov/cgi-bin/AtData/main_asd/ 



proved to be successful in separating accretion- 
powered vs. SF powered regions in the SINGS 
galaxy sample. 

We have calculated these ratios for our eight 
spectra, show n in Fig. [2] We have also plotted 
the data from Dale et al.1 (2006) that include nu- 
clear and extra-nuclear regions from SINGS, the 
SMC/LMC and galactic H II regions. Errors in the 
emission line fluxes calculation are around 10%. 
The boundaries in Fig. [21 are defin ed with the 
same curves as in iDale et al. (120061) . As we can 
see from Fig. El ULX-2, ULX-3, ULX-4, ULX- 

5, and ULX-6 are in zone I where Seyfert and 
Liners are observed, while ULX-1, comp-NW, and 
comp-CENT are located in zone III as star form- 
ing systems. For ULX-2 and ULX-6 the line ratio 
[S III]33.48 /.m/ [Si II]34.82 /im is < 0.01, which 
is extremely low compared with extragalactic ob- 
jects in Fig. El 

The SING S group has rece ntly developed mod- 
els (pahfit, ISmith et ai1l2007[) to characterize PAH 
features. The models include starlight contin- 
uum, featureless thermal dust continuum, pure 
rotational lines of H2, fine-structure lines, dust 
emission features, and dust extinction. The 
starlight continuum is represented by a black- 
body emission was fixed at = 5000 K. Nine 
thermal dust continuum components are rep- 
resented by blackbodies at fixed temperatures 
T m = 35,50,75,90,110,135,200,300,500 K. The 
program pahfit performs a fit through two hundred 
iterations, each iteration determines a x 2 value for 
all nine temperatures, until a minimum value for 
X 2 is found. Dust temperatures are based on 
the shape of the continuum spectrum. We show 
the corresponding fits in Table [TJ where we show 
ULX-2 and ULX-6 have hotter components, but 
the relative intensity of those components is low 
for ULX-1. 

The presence of a hotter dust component will 
flatten the spectra for our range of wavelengths. 
This effect is seen in Fig [T] for ULX-2 and ULX- 

6. The rest of the parameters modeled will be sub- 
ject of a forthcoming paper (Vazquez, et al. 2007b, 
in preparation). 

Finally, we have found a strong linear anti- 
correlation between the observed log(Lx) and 
log([S III]33.48 /mi/[Si II]34.82 /im). The regres- 
sion and intercept coefficients are —0.21 ±0.05 and 
39.09 ± 0.6, respectively with a RMS = 0.08. 
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4. Discussion 

For the first time the IR continuum and lines 
(~ 9.9 — 38 (im) have been observed in ULXs by 
using the power of the Spitzer IRS. It is likely 
that this ULXs are irradiating nearby star-forming 
clouds. Although the geometry is probably thus 
quite complex, we do find some hints of informa- 
tion on the underlaying X-ray radiation. The most 
striking result from the data analyzed in this work 
is the one in Fig. [2l which shows that five of the 
six ULXs reported bv lRoberts et al. ( 2002 ) fall in 
the "AGN zone" determined by the infrared diag- 
nostic emission line rat ios defined by th e emission 
lines of Ne, S, and Si jDale et al.ll2006l) . Both of 
the comparison regions used here, which do not 
harbor ULX sources, fall in the zone of SF. The 
SF zone is also occupied by extra-nuclear star- 
forming regions in the SINGS galaxy sample and 
SMC/LMC and galactic H II regions. The emis- 
sion line at [N I] 10.33 and other unidentified fea- 
tures also behave like [Si II] 34.82 //m, showing a 
stronger contribution in ULX-2 to 6 than in ULX- 
1 and both of the comparison regions. 

The presence of a hotter dust component in our 
PAH/continuum modeling (see Table [J), suggests 
that at least for ULX-2 and ULX-6 the [Si II] lines 
could be originating from an accreting source. A 
dilution factor could be affecting the rest of the 
ULXs possibly due to their location within their 
host galaxy. All of this is consistent with the anti- 
correlation we have found between the line emis- 
sion ratio of [S III]33.48/im/[Si II]34.82/xm and the 
observed X-ray luminosity of the six ULXs, the 
more luminous the sources the lower the ratio of 
these emission lines. 

Different scenarios for the strong [Si II] line 
have been suggested in the literature. In these 
scenarios the increasing density in H II, photo- 
dissociation (PDR), and X-ray dissociation re- 
gions (XDR), respectively, i ncreases the intensity 
of the [Si II]34.82 ixm line (seelMalonev. Hollenbach. 
1996[ Kaufman et"al1 l2006t iMeijerink k SpaansT 
20051 ). A third scenario suggests that heavy el- 
ements such as Si, Mg, and Fe may be returned 
to the gas phase by dust destruction (e.g., sput- 
tering) in regions subject to strong shocks caused 
by stellar winds, starbursts, and AGN activity 
( Dale et al. 20061 ). Details of the geometry and 
ionization balance in the region require more de- 



tailed fitting, and also spatial mapping of spectral 
parameters as ULXs are thought to have drifted 
from their birthplaces. The origin of these lines 
could thus be disentangled with detailed IRS map- 
ping observations and IRAC imaging, which are 
the subject of a forthcoming paper (Vazquez et 
al. 2007b). 

With the exception of a weak detection in 
ULX-1, we have not detected the emission line 
[O I V] 25.91 ^m, this line is a strong feature in SNR 



(e.g. iMorris et al.l2006l:IWilliams. Chu. & Gruendl 
20061 ) , supporting previous suggestions that ULX- 
1 is a S NR (e.g., based on radio and X-ray con- 
straints; iRoberts et al. I l2002l ). 

The mid-infrared emission line diagnostics used 
to uncover the nature of AGNs and SF galaxies 
has been shown here to work well separating SFRs 
and young SN from probably true a ccreting ULXs 



(acco rding to X-ray observations; IRoberts et al 
2002). But this is a small sample of ULXs to claim 



a general rule that could be improved observing a 
higher number of ULXs. The present work shows 
that a new diagnostic window has been opened on 
ULXs and their immediate environment via im- 
portant mid-infrared features like [Si II] 34.82 fim, 
made observationally available by the extraordi- 
nary performance of the IRS on board the Spitzer 
Space Telescope. 

Many thanks to Danny A. Dale who provided 
us his data for our figures. This work has been 
supported by the Spitzer grant 3360 (P.I. Colbert). 
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Fig. 1. — The combined spectra with the left SH 
and right LH modes for all of the sources analyzed 
in this study. For better visualization, the spec- 
tra have been plotted displaced along the y-axis. 
From ~ 9 to ~ 26.9 /xm spectra are displaced by 
0.75 Jy and ~ 26.9 to ~ 38 /im spectra are dis- 
placed by 15 Jy. Dashed lines are set to help in- 
dicate the slope of the continuum at wavelengths 
> 20 fim. 
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Fig. 2. — A neon, sulfur, and silicon diagnos- 
tic diagram showing line ratios for the eight re- 
gions in the NGC 4485/90 system (names in 
black). Nuclear and extra-nuclear reg ions from 
the sample of SINGS l|Dale et al.l I2006T) are also 
plotted (in red) together with SMC, LMC and 
galactic H II regions. The red line for ULX- 
2 means that the value for the line ratio of 
[Ne III] 15.55 /im/[Ne II] 12.81 /xm was not deter- 
mined due to the undetected lines of Ne for this 
source. Zones III and IV are typical of SF-regions. 
Objects i n zones I and II, c an be classified as AGN- 
powered ( Dale et al. 20061 ). 
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Table 1 

The list of ULXs of the NGC 4485/90 system taken from IRoberts et al.I (|2002h 



ULX 


RA (J2000) 


DEC (J2000) 


Notes 








T9. 7 










10 39 erg s" 1 


[K] 


[io- 10 ] 


[IO" 3 ] 


1 


12:30:29.5 


+41:39:27 


NGC 4490 Radio Src. Possible SNR 


1.0 (4.9) 


50, 75, 


72000, 13000 


3310 












200, 300 


0.1, 1.1 




2 


12:30:30.6 


+41:41:42 


NGC 4485 X-l 


4.0 (4.6) 


135, 200 


56, 2.3 


4810 


3 


12:30:30.8 


+41:39:11 


NGC 4490 Transient X-ray Src. 


2.9 (4.6) 


90, 110 


740, 240 


3.93 


4 


12:30:32.3 


+41:39:18 


NGC 4490 X-l 


1.9 (2.6) 


75, 90 


2500, 960 


1.75 


5 


12:30:36.3 


+41:38:37 


NGC 4490 X-2 


1.9 (2.6) 


75, 90 


8800, 16 


7.94 


6 


12:30:43.2 


+41:38:18 


NGC 4490 X-4 


3.1 (4.4) 


75, 90 


360, 580 


2990 












110, 135 


190, 10 





a Numbers in parenthesis give the intrinsic (unabsorbed) luminosity. 

b The total extinction as a function of the extinction at 9.7 /J,m (to. 7). The total extinction is modeled as a power low plus silicate 
features at 9.7 an d 18 fim. The features are counted as a combination of measured galactic extinction profiles for these wavelengths 
flSmith et aUI2007fl . 



